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Abstract 

We have successfully applied the Interacting Gluon Model (IGM) to calculate diffractive mass 
spectra measured recently in e-p collisions at HERA. We show that it is possible to treat them in 
terms of gluon-gluon collisions in the same way as was done before for hadronic collisions. Analysis 
of available data is performed. 

PACS number(s): 13.85.Qk, 11.55.Jy 



1 Introduction 

Diffractive scattering processes are related to the large rapidity gap physics usually interpreted in terms 
of Pomeron exchange Q . In hadronic diffractive scattering, one of the incoming hadrons emerges from 
the collision only slightly deflected and there is a large rapidity gap between it and the other final state 
particles resulting from the other excited hadron. In the standard Regge theory diffraction is visualised as 
due to the Pomeron exchange which implies that the excited mass spectrum behaves like 1 jM\ and does 
not depend on the energy g]. The same phenomena can, however, be understood without mentioning 

*e-mail: dunga@uspif.if.usp.br 
t e-mail: navarra@uspif.if.usp.br 
te-mail: wilk@fuw.edu.pl 



1 



the name of Pomeron by treating IP as a preformed colour singlet object consisting of only a part of the 
gluonic content of the diffractive projectile which is then absorbed by the other hadron ||. One possible 
realization of this idea was developed recently by us Q (for earlier similar attempts see ||). 

The first test of a theory (or a model) of diffractive dissociation (DD) is the ability to properly describe 
the mass (Mj) distribution of diffractive systems, which has been measured in many hadronic collision 
experiments Q and parametrized as (M x )~ a with a ~ f . 

Very recently diffractive mass spectra have been measured also in photoproduction processes at high 
energies at HERA Q . They were interpreted there in terms of standard Regge theory. In this work we 
would like to analyse these data using instead the Interacting Gluon Model (IGM) which was already 
successfully used to describe diffractive mass spectra and their energy dependences in hadronic reactions 
Q. The advantage of such approach is the possibility to check what part of the results is due to the 
simple implementation of conservation laws (notice that IGM was designed in such a way that the energy- 
momentum conservation is taken care of before all other dynamical aspects - a feature very appropriate 
for the study of all kinds of energy flows). It means that all notions to the Pomeron (or IP) in what 
follows is just symbolic representation of the DD process. 



2 IGM picture of a diffractive event 

As mentioned in Q], at the HERA electron-proton collider the bulk of the cross section corresponds to 
photoproduction, in which a beam electron is scattered through a very small angle and a quasi-real photon 
interacts with the proton. For such small virtualities the dominant interaction mechanism takes place 
via fluctuation of the photon into a hadronic state which interacts with the proton via the strong force. 
High energy photoproduction therefore exhibits similar characteristics to hadron-hadron interactions. 

In Fig. 1 we show schematically the IGM picture of a diffractive event in a photon-proton collision. 
According to it, during the interaction the photon is converted into a hadronic (mesonic) state and 
then interacts with the incoming proton ||. The meson-proton interaction follows then the usual IGM 
picture, namely: the valence quarks fly through essentially undisturbed whereas the gluonic clouds of 
both projectiles interact strongly with each other (by gluonic clouds we understand a sort of "effective 
gluons" which include also their fluctuations seen as qq sea pairs). The meson looses fraction x of its 
original momentum and gets excited forming what we call a leading jet (LJ) carrying xl = 1 — x fraction 
of the initial momentum. The proton, which we shall call here the diffracted proton, looses only a fraction 
y of its momentum but otherwise remains intact 

In the limit y — * 1, the whole available energy is stored in Mx which then remains at rest, i.e., Yx = 0. 
For small values of y we have small masses Mx located at large rapidities Yx- In order to regard our 
process as being trully of the DD type we must assume that all gluons from the target proton participating 
in the collision (i.e., those emitted from the lower vertex in Fig. 1) have to form a colour singlet. Only 
then a large rapidity gap will form separating the diffracted proton (in the lower part of our Fig. 1) and 
the Mx system (in its upper part), which is the experimental requirement defining a diffractive event. 
Otherwise a colour string would develop, connecting the diffracted proton and the diffractive cluster, and 
would eventually decay, filling the rapidity gap with produced secondaries. In this way we are effectively 
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introducing an object resembling closely to what is known as Pomeron (_ZP) and therefore in what follows 
we shall use this notion. The Pomeron may be treated as being composed of partons, i.e., gluons 
and sea qq pairs, in much the same way as hadrons, with some characteristic distribution functions which 
have been object of studies in HERA experiments JlJ, [T^ ]. 

As usual in the IGM Q we first start with the function x(x, y) describing the probability to form a 
central gluonic fireball (CF) carrying momentum fractions x and y of the two colliding projectiles: 

Xo 



• exp 



1 



[(y 2 )(x - (x)f + (x 2 )(y - (y)f - 2{xy)(x - (x))(y - (»))] 



(1) 



For our specific needs in this paper (application to DD events) where we are mostly interested in the 
x and M\ behaviour of the results, it is usefull to present ([!]) in the form where the x-dependence is 
factorized out: 
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Here \o denotes the normalization factor provided by the requirement that J^dx j^dy x(x,y)6(xy 

he 

possible CF and ^/s is proton-hadron center of mass energy. In the above expression y max = 



Kmin) = 1 with K m i n — ^| being the minimal inelasticity defined by the mass mo of the lightest 

This 

upper cut-off, not present in the non-diffractive formulation of the IGM (where y ma x = 1), is necessary 
to adapt the standard IGM to DD collisions. It is a kinematical restriction preventing the gluons coming 
from the diffracted proton (and forming our object IP) to carry more energy than what is released in 
the diffractive system. As it will be seen below, it plays a central role in the adaptation of the IGM to 
DD processes being responsable for its proper M\ dependence. The, so called, spectral function oj(x,y) 
contains all the dynamical inputs of the IGM in the general form given by (cf. ) 



w( ) = ^yf) G(x)G(2/) 0( xy . 

cr(s) 



Ki 



(5) 



where G's denote the effective number of gluons from the corresponding projectiles (approximated by 
the respective gluonic structure functions) and a gg and a are the gluonic and hadronic cross sections, 
respectively. In order to be more precise, the function G(y), as it can be seen in Fig. 1, represents 
the momentum distribution of the gluons belonging to the proton subset called Pomeron and y is the 
momentum fraction of the proton carried by one of these gluons. We shall therefore use the notation 
G(y) = Gp(y). This function should not be confused with the momentum distribution of the gluons 
inside the Pomeron, f g /p(P) (see below). 
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The moments (9™), q = x,y (we only require n = 1,2) are given by (Q) and are the only places where 
dynamical quantities like the gluonic and hadronic cross sections appear in the IGM. At this point we 
emphasize that we are all the time dealing with a meson (essentially the p°)-proton scattering. However, 
as was said above, we are in fact selecting a special class of events and therefore we must choose the correct 
dynamical inputs in the present situation, namely Gp(y) and the hadronic cross section a appearing in 
ui. 

As pointed out in the introduction, the Pomeron for us is just a collection of gluons which belong 
to the diffracted proton. In our previous work we have assumed that these gluons behave like all other 
ordinary gluons in the proton and have therefore the same momentum distribution. The only difference 
is the momentum sum rule, which for the gluons in IP is 

-1 

dyyG P {y) = p (6) 

where p ~ 0.05 (see [Q and below ) instead of p ~ 0.5, which holds for the entire gluon population in the 
proton. Alternatively we may treat the Pomeron structure in more detail and address the question of its 
"hardness" or "softness" . In order to make contact with the analysis performed by HERA experimental 
groups we consider two possible momentum distributions for the gluons inside IP : 

fit/r(P) = 6(1-/3) (7) 

and 

f' g ,M = q(J —/^ (8) 

where j3 is the momentum fraction of the Pomeron carried by the gluons and the superscripts h and s 
denote hard and soft respectively. We follow here the (standard) notation of ref. |l2). We shall use the 
Pomeron flux factor given by 

fp/ P {xp) = — (9) 

Xp 

where xp is the fraction of the proton momentum carried by the Pomeron and the normalization will 
be fixed later. Noticing that f3 = -p- the distribution Gp(y) needed in eq. (|) is then given by the 
convolution: 

g ' ( " ) - C-^'"""^' (i °> 

In our calculations we shall also use Gp{y) = 6 — ^ , the same expression already used by us before [Q. 
As it will be seen this choice corresponds to an intermediate between "soft" and "hard" Pomeron. In 
the upper leg of Fig. 1 we assume, for simplicity the vector meson to be p° and take G p ° (x) — G^(x). 
The fraction of diffracted nucleon momentum, p, allocated specifically to the P gluonic cluster and the 
hadronic cross section a are both unknown. However, they always appear in u as a ratio (— ) of parameters 
and different choices are possible. Just in order to make use of the present knowledge about the Pomeron, 
we shall choose 

a(s) = a Fp = a + b\n— (11) 

so 
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where So = 1 GeV 2 and a — 2.6mb and b = O.Olmb are parameters fixed from a previous systematic 
data analysis. As it can be seen, <j(s) turns out to be a very slowly varying function of yfs assuming 
values between 2.6 and 3.0 mb, which is a well accepted value for the Pomeron-proton cross section, 
and p ~ 0.05 (cf. Since the parameter £ has been fixed considering the proton-proton diffractive 

dissociation and we are now addressing the p — p° case we have some freedom to change a. In the 
following we shall also investigate the effect of small changes in the value of m on our final results. 

Although in the final numerical calculations the above complete formulation will be used, it is worth- 
while to present approximate analytical results in order to illustrate the main characteristic features of 
the IGM. It is straightforward to show that, keeping only the most singular terms in gluon distribution 
functions, i.e., G(x) — G h p{x) ~ \ and only the leading terms in y/s, one finds that for any choice of a gg 
in (|): 

(i) terms containing the mixed moment (xy) can be always neglected in comparison to those containing 
(x 2 ) or (y 2 ) (i.e., for example, D xy ~ (x 2 )(y 2 )); 



(ii) all moments are related in the following way: 

1 
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(x )~-(x); (y) ~ y max {x); (y ) ^ -y max {y)] 
i.e., all results can be expressed in terms of the (x) moment only; 
(Hi) the (x) moment has the following simple behaviour depending on the type of a gg chosen in (^J): 
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respectively. 

This allows us to write eq.(||) in a very simple form: 
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(14) 



term present in ( 14 ) can be traced back to the upper cut-off y = y n 



in (M) above. Because it is defined by the produced diffractive mass, y n 



A/T 



it provides then 



automatically the -^y behaviour. The other two factors have a much weaker dependence on M\ and 
they tend to compensate each other (they provide, however, all possible non-trivial energy dependence 
ofDD, cf. |). 



3 Comparison with experimental data 

The IGM diffractive mass spectrum is given by Q : 

dN r 1 r 1 



, w . , dx I dy x(x, y) 5 (M|- — sy) 9 (xy — K 2 nm ; 

uivi x Jq Jq 

-! dx X (x,y=^pj, (15) 



5 



or, in the approximate form, 



6M\ Mf n{x) 



exp 



, dx exp 



(16) 



We would like to emphasize two aspects of the approximate formula above: 



(i) it explicitely exhibits the characteristic M\ dependence of diffractive collisions, namely it shows 
the [M\)~ a behaviour with a ~ 1; 

(ii) the exponential and integral factors have a very weak dependence on M\ but they contain a non- 
trivial energy yfs dependence, which is intrinsic to the model and comes ultimately from phase 
space limits and cross sections contained in eq.(^|). 

In Fig. 2 we compare eq.([l5|) with the recent data from the HI collaboration [Q. In all formulas yfs 
will now be replaced by W, the photon-proton center of mass energy. Figure 2a (2b) presents data for 
W = 187 GeV (W = 231GeV). The different curves correspond to the choices I (mo = 0.31 GeV, a = 2.7 
mb), II (m = 0.35 GeV, a = 2.7 mb), III (m = 0.31 GeV, a = 5.4 mb) and IV (m = 0.35 GeV, a = 5.4 
nib), respectively. In all these curves we have used Gp{y) — 6^—^-. As expected, the distribution at 
low M\ is very sensitive to threshold effects. When we go from the upper to the lower solid (dashed) 
lines we can observe that the increase of the Pomeron-hadron cross section changes the distribution in 
such a way that larger masses M\ are favoured. 

We would like to stress that in curve II there is no free or new parameter. All parameter values are 
the same as in our previous paper devoted to hadronic diffraction. It misses only the very small mass 
region points, where we expect it to be below the data, since we do not include resonance effects. In the 
large mass region a better agreement with data may be achieved with a somewhat larger value of the 
Pomeron-hadron cross section. This region may, however, be influenced by other effects, one of which we 
discuss below. 

In Fig. 3 we compare the same data (W — 187 GeV) with our mass spectrum obtained with G P {y) 
(curve I), Gp(y) (curve II) and G P {y) (curve III). This comparison suggests that the "hard" Pomeron 
can give a good description of data. The same can be said about our "mixed" Pomeron, which, in fact 
seems to be more hard than soft. These three curves were calculated with exactly the same parameters 
and normalizations, the only difference being the Pomeron profile. Apparently the "soft" Pomeron (curve 
III) is ruled out by data. Curve IV shows, however, that with a different choice of parameters mo = 0.50 
GeV and a = 5.4 mb a good agreement is again obtained. Considering the large ammount of data already 
described previously by the IGM, this choice is extreme. We conclude therefore that the "soft" Pomeron 
is disfavoured. This same conclusion was found in refs. |l2|, |l3[ . 

A very interesting question regarding DD processes is whether or not semihard interactions play a 
role in diffractive physics. In hadronic non-diffractive collisions semihard scatterings are expected to 
be visible at c.m.s. energies around yfs ~ 500 GeV. In such scatterings two partons interact with a 
momentum transfer of pr — 2 — 4 GeV, forming two so-called minijets. Since Kqcd "C Pt -C v^i 
minijet cross sections can be calculated with perturbative QCD and they are large enough to be relevant 
for minimum bias physics. In the IGM, energy deposition is occuring due to gluon-gluon collisions in 
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both perturbative (semihard) and non-perturbative regimes. The gluon-gluon cross section in eq. (|||) is 
computed with perturbative QCD or with a non-perturbative ansatz according to the scale (= xys). The 
relative importance of minijets with respect to the soft processes was fixed following the experimental 
estimates of the minijet cross section made by the CERN UAl and UA5 collaborations in hadronic 
collisions. We have been assuming so far that the onset of semihard physics in DD occurs at the same 
energy as in ordinary non-diffractive processes. This may not be true, or even if it is true, there are 
uncertainties regarding the precise value of the relevant energy scale. In Fig. 4 we repeat the fit of 
Fig. 2 using only curves II, which are our "conservative" calculation, plotted with solid lines. Since 
the present energies are not yet very large, we could just neglect the small minijet component. The 
dashed curves in Fig. 4 show the effect of switching off the semihard contribution. There is only a small 
enhancement in the tail of the spectra. Without minijets the energy deposition in the central blob of Fig. 
1 is decreased and the leading particle, in the upper leg (p meson after interaction), is more energetic. It 
contributes more to the diffractive mass M\ and makes it larger. This effect is negligible for very low M\ 
but becomes visible at large diffractive masses. Repeating this comparison (total spectrum versus the 
spectrum without minijets) at higher energies we observe that the magnitude of the minijet contribution 
is always small and shows always the tendency to produce slightly faster falling distributions at the end 
of the spectrum. This suggests that in DD processes minijets are unimportant even at very high energies. 



4 Summary and conclusions 

In conclusion, a straightforward (and with no new parameter) extension of our model of hadronic diffrac- 
tion to photon-proton reactions is able to fit the data for diffractive mass excitation presented in ref. Jt| 
within small discrepancies. The agreement may become better with some small changes motivated by 
uncertainties in previous fitting procedures. 

Our analysis of data suggests that the Pomeron is "hard" and undergoes mostly "soft" interactions. 
This means that this object is composed by a relatively "small" number of gluons carrying each, in the 
average a large fraction of the Pomeron momentum, but a small fraction of the total momentum of the 
proton and undergoing mostly soft (with respect to a hard energy scale ~ 2 — 3 GeV) collisions with the 
gluons of the other hadron. 

The fact that our model is successful means that energy flow in many and different high energy 
reactions can be understood as an incoherent superposition of parton-parton scatterings constrained by 
energy conservation. 
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Figure Captions 

Fig. 1 IGM description of a photon-proton scattering with the formation of a diffractive system of invariant 
mass Mx- 

Fig. 2a Diffractive mass spectrum for jp collisions at W = 187 GeV calculted with the IGM (eq.(|l5|)) and 
compared with HI data 0. The different curves correspond to the choices: I (mo = 0.31 GeV, 
a = 2.7 mb), II (m = 0.35 GeV, a = 2.7 mb), III (m = 0.31 GeV, a = 5.4 mb) and IV (m = 0.35 
GeV, a = 5.4 mb), respectively. 

Fig. 2b The same as Fig. 2a for W = 231 GeV. 

Fig. 3 Data from ref. compared with eq. jlif). The solid line (curve II) corresponds to the choice 
mo = 0.35 GeV, a — 2.7 mb and Gp(y). Curves I (dashed) and III (dotted) are obtained replacing 
Gp(y) by G P (y) and G s P (y) respectively. Curve IV is obtained with G s P (y) and m = 0.50 GeV 
and a = 5.4 mb. 
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Fig. 4a Diffractive mass spectrum for collisions at W = 187 GcV, curve II of Fig. 2a, shown with a solid 
line and compared with the same spectrum without the minijet contribution (dashed line). 

Fig. 4b The same as Fig. 3a for W = 231 GeV. 
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Figure 1 
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Figure 3 
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